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INTRODUCTION 
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This  report  presents  the  results  of  a study  of  detectors  for 
long-period  bodywave  phases.  It  was  motivated  by  the  fact  that  there  is  a 
good  separation  between  the  earthquake  and  explosion  populations  by  means 
of  long-period  bodywave  magnitudes  (Kisslinger  et  al.  , 1974),  and  these 
phases  therefore  provide  additional  discriminatory  power  where  they  can  be 
detected. 

When  a simple  bandpass  filter  is  used  as  the  only  processing 
scheme  the  50%  detection  threshold,  the  short-period  magnitude  at  which  half 
of  all  events  are  detected  by  means  of  long-period  bodywaves,  is  about  m^  = 

5.  5 for  P waves  and  about  5.  1 for  S waves,  considerably  higher  than  the  cor- 
responding threshold  for  surface  waves  of  about  4.  6 (Strauss,  1977). 

The  processors  described  in  this  report  take  advantage  of  the 
polarization  of  P and  S waves.  Test  events  show  good  improvement  in  signal- 
to-noise  ratio,  but  when  evaluated  on  a data  set  consisting  of  earthquakes  in 
the  Kurile  Islands  recorded  at  the  Guam  Seismic  Research  Observatory,  no 
reduction  in  threshold  was  achieved.  There  is  evidence  that  this  is  due  to  the 
particular  data  set  used,  and  that  the  P detection  threshold  may  be  as  low  as 
m = 5. 0. 

D 

In  Section  II  the  processing  schemes  used  here  are  described, 
and  examples  of  their  effects  on  data  presented.  Section  HI  describes  the  data 
sample  investigated,  and  Section  IV  gives  detection  thresholds  and  long-period 
bodywave  magnitudes  found  here.  Section  V presents  conclusions  and  recom- 
mendations. 
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SECTION  n 
PROCESSORS 

A natural  first  step  in  the  extraction  of  long-period  bodywaves 
from  seismic  noise  is  to  form  beams  at  the  appropriate  velocity  at  arrays 
such  as  the  Iranian  Long  Period  Array  (ILPA)  and  the  reconfigured  ALPA  ar- 
ray in  Alaska.  This  approach  was  discarded  for  the  following  reasons. 

Since  ILPA  and  the  reconfigured  ALPA  are  only  seven- element 
arrays,  the  maximum  noise  suppression  to  be  found  there  assuming  uncor- 
related noise  is  about  8 dB.  Strauss  (1976)  has  found  signal-to-noise  ratio 
improvements  on  the  order  of  3 to  6 dB  there,  consistent  with  the  expected 
suppression  of  uncorrelated  noise  and  some  signal  degradation.  His  results 
apply  to  surface  waves  propagating  at  3.  6 km/sec.  In  the  present  study  beams 
on  a few  large  events  at  the  appropriate  bodywave  phase  velocity  of  12.  5 km/ sec 
did  not  yield  even  those  reduced  gains.  This  may  be  attributed  to  two  causes. 

First,  visual  examination  of  individual  site  records  indicates 
that  site-to-site  similarity  is  not  very  good  for  long-period  bodywaves.  Such 
records  are  presented  in  Figure  II- 1 for  an  event  from  the  Greece- Albanian 
border.  Considerable  signal  amplitude  degradation  may  be  expected  in  beams 
from  such  an  event. 

Second,  the  lack  of  similarity  displayed  in  Figure  II- 1 suggests 
that  beams  at  P and  S velocities  will  not  only  degrade  amplitudes  but  wave- 
forms, in  the  sense  that  the  plane  and  linear  polarization  of  these  phases  will 
be  contaminated  by  beamforming.  Since  previous  experience  with  signal  pro- 
cessors has  suggested  that  adaptive  filters  searching  for  constants  of  the 
motion  are  highly  effective  (Lane,  1976),  it  was  decided  to  use  data  in  which 
signals  were  best  preserved. 
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Long-period  data  recorded  at.  the  Seismic  Research  Observa- 
tories (SRO)  were  chosen  for  this  purpose.  These  results  will  then  be  ap- 
plicable to  this  increasingly  important  type  of  data,  and  can  be  extrapolated 
to  beam  data  if  beams  which  preserve  signal  similarity  can  be  found. 

Since  adaptive  filters  are  non-linear  due  to  their  dependence 
on  the  data,  it  is  important  that  they  be  preceded  by  the  best  possible  prefilter. 
This  will  be  in  general  a Wiener  filter,  using  as  a reference  spectrum  one 
from  a typical  event  from  the  region  of  interest,  or  an  average  over  such 
spectra.  However,  examination  of  the  few  large  events  from  the  Kurile 
Islands,  the  region  used  in  this  study,  showed  that  bodywaves  from  there  show 
as  much  variation  between  one  another  as  they  do  from  site  to  site  as  shown 
in  Figure  11- 1.  An  average  of  such  signals  would  probably  not  be  representa- 
tive of  the  Kurile  Island  earthquake  population.  Consequently,  a simple  band- 
pass filter,  rejecting  all  energy  outside  the  frequency  range  0.024-0.  059  Hz, 
was  used  as  the  prefilter  in  this  study. 

The  adaptive  filters  used  here  searched  for  linearly  polarized 
motion  in  the  correct  emergence  direction  for  extracting  P waves  and  plane 
polarized  motion  perpendicular  to  the  emergence  direction  to  find  S waves. 
Following  the  form  of  the  three  component  adaptive  filter  developed  for  surface 
waves  (Lane,  1976),  each  filter  segments  the  vertical,  radial,  and  transverse 
motion,  Fourier  transforms  each  segment,  applies  a filter  whose  weights  in 
the  frequency  domain  depended  on  the  behavior  of  the  total  motion,  inverse 
transforms,  and  smooths  together  overlapping  segments  to  produce  a continu- 
ous output  waveform. 

Variations  in  the  design  of  such  filters  lie  in  the  way  the  filter 
weights  depend  on  the  data.  In  the  present  case  two  approaches  were  tried. 
First,  each  Fourier  coefficient  within  the  frequency  interval  where  the  signal 
energy  was  expected  was  tested  separately  for  the  correct  polarization,  and 
an  independent  weight  was  assigned  to  real  and  imaginary  coefficients. 
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Because  such  a filter  is  not  zero  phase,  and  disperses  the  signal,  it  was  re- 
jected during  preliminary  testing  in  favor  of  an  approach  where  the  weights 
were  determined  from  and  applied  to  the  total  Fourier  amplitude. 


To  sear  h for  linearly  polarized  motion,  a spatial  propagation 
vector  P was  formed  at  each  frequency  as  shown  in  Figure  11-2.  Its  ampli- 
tude along  each  axis  was  the  square  root  of  the  sum  of  the  squares  of  imagin- 
ary and  real  Fourier  coefficients  of  the  motion  along  that  axis,  suitably  nor- 
malized so  that  the  vector's  amplitude  was  unity.  The  scalar  product  between 
this  vector  and  a unit  vector  ^ in  the  expected  propagation  direction,  as  de- 
termined by  the  epicentral  distance,  equals  the  cosine  of  the  angle  9 between 
the  apparent  propagation  direction  of  energy  at  that  frequency  and  the  propaga- 
tion direction  for  P waves  from  the  region  of  interest.  In  general  P will  not 
lie  in  the  vertical- radial  plane,  due  to  contamination  of  the  signal  with  noise. 


If  0 is  small,  the  energy  may  be  attributed  to  a P wave,  where- 
as if  it  is  near  90°,  it  might  be  an  S wave  which  have  nearly  identical  emer- 
gence directions  to  P waves.  Filter  weights  for  extraction  of  either  P or  S 
waves  can  then  be  assigned  using  the  same  angle  as  a criterion.  Choosing 
the  way  in  which  to  do  this  is  the  remaining  part  of  the  filter  design. 


A Monte  Carlo  method  was  used  to  find  the  probability  distribu- 
tion function  of  the  angle  between  the  true  and  apparent  propagation  directions, 
assuming  a model  with  a P wave  in  the  true  propagation  direction,  and  addi- 
tive uncorrelated  noise.  This  procedure  has  been  discussed  more  fully  in 
Lane  (1976).  The  probability  distribution  found  for  narrowband  signal-to- 
noise  ratios  near  2.  0 may  be  roughly  approximated  by  the  curve  of  Figure  II- 3, 
which  was  chosen  for  the  filter  weighting  function.  At  higher  signal- to-noise 
ratios  the  probability  distribution  function  peaks  more  sharply  near  zero  de- 
grees, while  at  low  signal-to-noise  ratios  is  becomes  flat.  The  probability 
distribution  function  is  of  course  normalized  to  unit  area,  while  the  filter 
weights  of  Figure  II- 3 are  chosen  to  be  one  at  angles  less  than  20  degrees. 
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The  same  filter  weights  were  chosen  for  the  S wave  filter,  ex- 
cept that  in  that  case  the  angle  plotted  on  the  ordinate  of  Figure  II- 3 refers 
to  the  angular  distance  of  the  particle  motion  vector  from  a plane  perpendic- 
ular to  the  propagation  vector,  i.  e.  , ff/2-0  . Thus,  both  SH  and  SV  waves 
will  be  passed  by  this  filter. 

Only  the  vertical  and  radial  traces  are  displayed  for  the  output 
of  the  P wave  filter,  since  no  motion  is  expected  on  the  transverse  component. 
All  three  components  of  the  S wave  filter  are  displayed,  since  in  general  both 
SV  and  SH  waves  are  possible.  Such  an  output  for  a large  event,  for  which  P 
and  S phases  were  detected  on  the  bandpass  filtered  traces,  is  shown  in  Figure 
II-4.  The  P wave  filter  traces  do  not  extend  past  the  S wave  arrival  time, 
since  S motion  often  passed  through  the  P wave  filter  and  dominated  the  trace. 
The  complementary  effect  can  be  seen  on  the  S wave  traces  in  Figure  II-4  be- 
tween the  P and  S phase  arrival  times,  where  there  are  large  deflections  pre- 
sumably due  to  P motion  leaking  through  the  S filter.  A smaller  event  is 
shown  in  Figure  II- 5,  and  was  detected  on  both  P and  S filters.  Considerable 
noise  precedes  this  event. 

Leakage  of  compressional  wave  energy  through  the  sheer  wave 
filter,  to  appear  as  motion  between  the  P and  S arrival  times  on  the  S wave 
detector,  implies  that  the  models  used  here  do  not  fit  the  data  particularly 
well.  This  is  tirue  for  both  large  and  small  events,  where  it  is  a more  serious 
problem,  since  it  suggests  that  the  design  of  each  detector  could  be  improved. 
Such  improvement  should  be  part  of  any  further  study  of  this  problem. 

Noise  suppression  for  the  event  in  Figure  II- 4 was  complete, 
in  that  the  RMS  noise  amplitude  before  the  signal  arrival  was  zero.  The  noise 
suppression  by  the  S wave  filter  ranged  between  25  and  35  dB.  For  the  event 
in  Figure  II-5,  noise  suppression  was  much  less,  as  were  starting  noise  am- 
plitudes. The  P wave  filter  suppressed  the  noise  by  10  to  15  dB,  while  the 
suppression  by  the  S wave  filter  ranged  from  5 to  15  dB.  These  figures  imply 
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gains  from  0.  2 to  0.  7 magnitude  units  if  signals  are  passed  unchanged.  Signal 
degradation  amounts  to  about  5 dB,  however,  presvimably  caused  again  by  lack 
of  agreement  between  model  and  data. 
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SECTION  III 
DATA  SET 

Since  the  most  convincing  test  of  a detector  is  to  evaluate  its 
performance  directly  on  a large  suite  of  events,  that  procedure  was  used  in 
this  study.  All  earthquakes  reported  by  NORSAR  between  1 January  1976  and 
27  June  1976  in  the  Kurile  Islands  region  were  examined  with  the  detector 
described  in  Section  II.  NORSAR  were  used  throughout  this  study.  It 

was  hoped  that  the  combination  of  constant  event  location  and  recording  station 
would  remove  the  effect  of  those  variables  from  the  evaluation  of  the  detector's 
reduction  in  detection  threshold. 

Guam  is  32°  distant  from  the  Kurile  Islands,  which  form  an 
active  seismic  region.  This  distance  is  a reasonable  one  at  which  to  search 
for  bodywaves  in  that  the  P and  S phases  are  well  separated  in  arrival  times, 
and  yet  may  still  have  reasonably  large  amplitudes.  A total  of  114  events 
were  available,  ranging  in  short-period  bodywave  magnitude  from  3.  3 to  5.  7. 

Some  large  events  in  this  set  were  found  to  be  part  of  an  earth- 
quake swarm,  and  consequently  the  root-mean- square  noise  levels  at  the  pre- 
dicted P and  S wave  arrival  times  for  most  of  these  events  were  as  high  as 
800  lim,  roughly  50  times  the  average  value  when  no  swarm  was  taking  place. 
Since  the  model  developed  by  Ringdal  (1974)  and  used  here  for  estimating  de- 
tection thresholds  assumes  a Gaussian  distribution  of  noise  magnitudes,  those 
events  arriving  in  a background  of  systematically  rather  than  randomly  vary- 
ing noise  cannot  be  legitimately  used  for  threshold  estimation.  Accordingly, 
events  with  noise  amplitudes  preceding  the  signal  arrival  of  more  than  100 
Mm  were  eliminated  from  the  data  set,  leaving  a total  of  100  events.  Their 
distribution  is  shown  in  Figure  III-l. 
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SECTION  IV 
RESULTS 


All  100  events  were  examined  on  the  bandpass  filtered  traces 
for  P or  S detections.  Detection  was  allowed  if  the  peak  motion  occurred 
within  40  seconds  of  the  predicted  arrival  time,  and  that  peak  was  twice  as 
large  as  any  peak  in  the  preceding  600  seconds  of  the  record.  Detections 
claimed  for  bodywaves  will  not  be  as  reliable  as  those  for  surface  waves, 
however,  due  to  the  lack  of  dispersion  in  bodywaves,  an  important  require- 
ment for  detection  of  surface  waves. 

A.  DETECTION  THRESHOLDS 

Detection  histograms  and  best  fits  of  Ringdal's  (1975)  constant 
plus  cumulative  Gaussian  distribution  to  the  associated  detection  probabilities 
are  shown  for  the  bandpass  filtered  traces  in  Figures  IV-1  and  IV-2  for  P and 
S phases,  respectively.  The  data  for  P waves  is  so  sparse  that  it  can  only  be 
claimed  that  the  50%  detection  threshold  is  high.  The  data  fit  the  model  much 
better  for  S detections,  where  the  50%  detection  threshold  is  5.  1. 

Next,  all  100  events  were  processed  with  the  standard  filter 
weights,  and  detections  claimed  using  the  criteria  listed  above,  with  one 
modification.  The  interval  between  the  P and  S arrival  times  on  the  S wave 
traces  were  disregarded  when  requiring  that  the  S motion  be  twice  the  ampli- 
tude of  a peak  in  the  preceding  noise  gate,  and  the  same  noise  gate  was  used 
for  P and  S detections.  This  eliminated  the  possibility  that  valid  S detections 
might  be  rejected  due  to  leakage  of  P waves  through  the  S detector  noise  gate, 
as  occurred  for  the  large  event  of  Figure  II- 3. 
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Detection  histograms  and  the  maximum  likelihood  fit  to  the  de- 
tection probabilities  are  shown  for  P and  S waves  in  Figures  IV- 3 and  IV- 4. 

Despite  the  increased  number  of  detections  at  low  magnitudes  over  the  band- 
pass filtered  results,  the  detection  threshold  has  been  increased  at  the  90% 
detection  level  for  both  phases.  The  distribution  of  detections  for  the  P de- 
tector is  such  that  they  are  interpreted  as  false  alarms,  while  one  event  has 
been  lost  from  detected  status.  The  S detector  produced  the  same  number  of 
detections,  but  at  different  magnitudes,  as  did  the  bandpass  filter,  raising  its 
threshold  slightly. 

In  an  effort  to  improve  these  results,  all  100  events  were  re- 
run using  different  filter  weighting  functions.  First  a function  which  was  zero 
if  its  argument  was  greater  than  18°  and  one  if  less  was  tested,  corresponding 
to  high  signal-to-noise  ratio  events.  Another,  which  searched  for  low  signal- 
to-noise  ratio  events,  rejected  motion  more  than  35°  from  the  predicted  di- 
rection and  accepted  all  others.  In  both  cases  the  results,  in  terms  of  detec- 
tion thresholds,  were  markedly  inferior  to  the  standard  filter  weights  derived 
from  the  Monte  Carlo  technique. 

Finally,  the  results  of  processing  with  the  standard  filter 
weights  were  reexamined  using  much  less  strict  detection  criteria.  If  there 
was  a distinct  pulse  of  energy  near  the  predicted  arrival  time  which  looked 
'seismic',  i.  e. , consisted  of  only  a few  cycles,  a detection  was  allowed.  This 
procedure  lead  to  a high  false  alarm  rate,  but  this  was  allowed  for  in  the 
maximum  likelihood  fitting  process.  Figures  IV- 5 and  IV- 6 show  cumulative 
Gaussian  fits  to  the  detection  probabilities  with  a variable  false  alarm  rate, 
also  adjusted  to  maximize  the  likelihood  function. 

The  50%  and  90%  detection  thresholds  listed  in  the  lower  right 
hand  corner  of  these  plots  are  the  correct  thresholds,  rather  than  the  point  i 

at  which  the  detection  probability  reaches  0.  5 and  0.  9.  The  latter  figures 
include  the  contribution  of  false  alarms,  whereas  the  former  are  parameters 
in  the  fit  to  true  detections. 
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A substantial  reduction  in  the  P detection  threshold  has  been 

i accomplished  at  the  expense  of  a false  alarm  rate  of  about  0.  26.  The  50% 

• « 

threshold  at  this  alarm  rate  is  about  m = 4.  8,  down  from  5.  5 without  the  de- 

b 

I tector.  The  S detector  shows  no  real  improvement  over  the  bandpass  filter, 

even  though  its  false  alarm  rate  increased  to  0.  15. 

I In  Figure  IV-7  the  epicenters  of  all  those  events  in  the  data 

sample  for  which  PDE  data  were  available  are  plotted.  The  concentration  of 
[ , events  near  44.  5°N  and  149.  4°E  suggests  that  those  events  share  a common 

mechanism.  That  mechanism  may  be  one  which  results  in  low  bodywave 
; . magnitude  for  P waves  at  Guam,  since  all  but  two  of  these  events  were  unde- 

tected with  the  P detector.  Furthermore,  those  two  events  of  m = 5,  7 which 

b 

i were  detected  had  long-period  P wave  magnitudes  a full  magnitude  unit  lower 

than  that  of  the  third  5.  7 event  in  the  sample,  which  was  also  detected 
but  lay  outside  the  region  in  question  at  the  point  marked  'A'  in  Figure  IV-7. 

Reliance  on  NORSAR  detections  to  develop  the  data  set  has 
probably  aggravated  this  problem.  It  is  plausible  that  all  the  events  in  the 
data  set  might  have  lobes  in  their  P wave  radiation  patterns  toward  NORSAR 
but  nulls  near  the  direction  of  Guam.  In  that  case  no  detector  would  give 
good  results.  A much  more  satisfactory  data  set  would  be  one  which  depend- 
ed for  detections  and  short-period  magnitudes  on  the  PDE  list  or  on  short- 
period  data  taken  at  the  same  site  as  the  long-period  data. 

B.  MAGNITUDES 

Despite  these  discouraging  results,  long-period  bodywave 
magnitudes  were  calculated  for  all  events  detected  using  the  criteria  which 
led  to  the  histogram  of  Figures  IV-3  and  lV-4.  The  standard  expression  for 
short-period  bodywave  magnitudes  was  used, 

% = ("t  ) ^ ^ 
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where  A is  the  peak  amplitude  of  the  vertical  trace  in  millimicrons,  T the 
period  in  seconds,  and  C a constant  dependent  on  distance  and  tabulated  by 
Vieth  and  Clawson  (1972).  These  magnitudes  are  plotted  versus  short-period 


magnitude  in  Figure  IV- 8 for  P detections. 

If  we  accept  the  hypothesis  that  the  events  in  the  clustered  near 

44.  5°N  and  149.  5°E  are  part  of  a swarm  with  low  long-period  amplitudes  at 

Guam  and  eliminate  them  from  Figure  IV- 8,  only  the  two  points  at  2.  9, 

m = 5.  7 are  dropped.  The  solid  line  is  the  best  linear  fit  to  the  rest  of  the 
b 

data.  Considering  the  large  number  of  false  alarms  implied  by  the  detection 
curve  of  Figure  IV- 3,  the  correlation  coefficient  of  0.  6 is  probably  reasonable. 
Since  corner  frequencies  for  earthquakes  in  this  magnitude  range  are  greater 
than  the  frequency  at  which  short-period  bodywave  magnitudes  are  measured 
(Sax,  1975),  both  long-period  and  short-period  magnitudes  are  measured  on 
the  flat  part  of  the  earthquake  source  spectrum.  Consequently,  the  slope  of 
the  best  fit  to  the  data  in  Figure  IV-8  should  be  1.  0,  and  the  slope  found  there 
of  0.  64  is  not  unreasonable  considering  the  low  quality  of  the  data.  Because 
of  the  large  variance  apparent  in  the  m^  measurements  of  P waves,  consider- 
able more  data  are  needed  to  establish  their  validity  as  discriminant  measure- 
ments. 

In  contrast  to  the  P detector  results,  all  but  one  large  event 
from  those  detected  by  the  S wave  detector  lie  within  the  earthquake  swarm 
region  near  44.  5°N,  149.  5°E,  as  shown  in  Figure  IV-9.  This  is  not  unreason- 
able, since  P and  S waves  do  not  have  the  same  radiation  pattern.  Further- 
more, all  three  magnitude  5.  7 events  have  about  the  same  long-period  body- 
wave  magnitude,  and  this  magnitude  is  consistent  with  that  of  the  smaller 
events,  as  shown  in  Figure  IV-IO,  where  m^^  is  plotted  versus  for  S waves. 
The  slope  of  the  best  fit  to  the  data  is  0.  53,  consistent  with  that  found  for 
long-period  P waves.  Thus,  the  S wave  data  as  well  as  the  P wave  data  are 
consistent  with  the  assumption  that  the  particular  data  sample  used  here  con- 
tributed to  the  poor  results  obtained.  If  this  interpretation  is  correct,  it  is 
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important  to  evaluate  any  future  detector  design  on  a data  sample  with  wide 
enough  spatial  distribution  that  no  bias  due  to  radiation  patterns  is  introduced. 


The  rather  large  values  of  sigma  found  on  some  of  the  fits  to 
detection  probabilities  are  probably  due  in  part  to  this  separation  of  the  data 
into  two  populations.  This  results  in  wide  variations  in  the  long-period  body- 
wave  magnitude  at  a given  short-period  magnitude,  and  sigma  is  a measure- 
ment, in  part,  of  this  variation. 

We  may  conclude  from  these  results  that  the  P detector  at 
least  is  behaving  qualitatively  as  it  should,  in  that  it  passes  signals  while 
rejecting  noise,  but  that  the  gain  so  achieved  is  not  sufficient  to  make  clear 
detections  possible.  The  S detector  seems  to  be  behaving  much  worse,  prob- 
ably due  to  poor  agreement  between  data  and  model,  as  evidenced  by  leakage 
of  S energy  onto  the  P detector  traces. 


SECTION  V 
CONCLUSIONS 


Two  conclusions  can  be  drawn  from  this  study.  First,  while 
the  detectors  studied  here  appear  to  be  promising  in  terms  of  improvement 
in  detection  capability  of  long-period  P waves,  leakage  of  energy  from  phases 
which  should  be  nulled  by  the  detector  probably  leads  to  less  than  optimum 
performance.  A model  for  the  motion  of  P waves  closer  to  the  true  particle 
motion  should  reduce  this  leakage  and  lower  the  detection  threshold.  Also, 

S wave  particle  motion  should  be  investigated  to  explain  the  lack  of  positive 
results. 

Second,  there  is  evidence  that  part  of  the  responsibility  for  the 
P detectors  performance  lies  in  apparent  nulls  in  the  radiation  pattern  peculiar 
to  the  data  set  used  here.  If  this  is  the  case,  another  data  set  should  be  used 
for  any  future  evaluations.  Some  means  may  also  be  needed  to  correctly  ac- 
count for  the  effects  of  multiple  radiation  patterns  on  the  detection  capability 
of  events  from  a given  source  region.  Some  consideration  should  be  given  to 
a joint  measurement  of  long-period  P wave  and  S wave  bodywave  magnitude 
which  minimizes  the  effect  of  radiation  patterns  on  detection  and  measurement 
of  magnitude. 
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